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Distorted black holes radiate gravitational waves. In the so-called ringdown phase radiation is 
emitted as a discrete set of complex quasinormal frequencies, whose values depend only on the black 
hole's mass and angular momentum. Ringdown radiation could be detectable with large signal-to- 
noise ratio by the Laser Interferometer Space Antenna (LISA). If more than one mode is detected, 
tests of the black hole nature of the source become possible. The detectability of different modes 
depends on their relative excitation, which in turn depends on the cause of the perturbation (i.e. on 
the initial data). A "universal", initial data-independent measure of the relative mode excitation is 
encoded in the poles of the Green's function that propagates small perturbations of the geometry 
("excitation factors"). We compute the excitation factors for general-spin perturbations of Kerr 
black holes. We find that for corotating modes with I — m the excitation factors tend to zero in the 
extremal limit, and that the contribution of the overtones should be more significant when the black 
hole is fast rotating. We also present the first analytical calculation of the large-damping asymptotics 
of the excitation factors for static black holes, including the Schwarzschild and Reissner-Nordstrom 
metrics. This is an important step to determine the convergence properties of the quasinormal mode 
expansion. 

PACS numbers: 04.70.-s, 04.30.Db, 04.80.Cc, 04.80.Nn 



I. INTRODUCTION 



Distorted black holes emit gravitational radiation as a discrete sum of quasinormal modes (QNMs), damped os- 
cillations whose frequencies and damping times depend only on the black hole's mass and angular momentum. The 
QNM frequencies scale with the inverse of the black hole's mass. Therefore, the optimal sensitivity bandwidth of 
each detector determines the mass range that can be detected. The dominant QNM frequency for quadrupolar ra- 
diation from a Schwarzschild black hole is / = 1.207 x 10-2(10^Mq/M) Hz. As a consequence, the collapse of 
Population-Ill stars with mass M > 26OM0 or 25Mq < M < AOMq forms massive black hole remnants that could be 
detectable by ground-based, high-frequency gravitational wave interferometers such as the advanced Laser Interferom- 
eter Gravitational- Wave Observatory (LIGO) [l| and Virgo [^]. The space-based Laser Interferometer Space Antenna 
(LISA), being sensitive in the frequency band ~ 10^^ — 10~^ Hz, will detect the gravitational radiation emitted by 
oscillating black holes of mass IO^Mq < M < 10^ — 10^ Mq with large signal-to-noise ratio (SNR) throughout the 
observable Universe [3, 3- With large SNR come precise measurements of the source parameters, and an extraordinary 
opportunity to study the physics of massive black holes Q. 

Since the LISA SNR can be very large, more than one mode could be detected in the ringdown waveform. Such a 
detection would allow an unprecedented test of the black hole nature of the source. The basic idea is quite simple. 
In general relativity the complex QNM frequencies form a discrete set uJimn classified by three integers: the indices 
{I, m) come from the separation of the angular dependence of the perturbations, and the index n labels frequencies by 
the magnitude of their imaginary part (large n means large imaginary part and short damping time). Because of the 
"no-hair theorem" , QNM frequencies depend only on the mass and angular momentum of a Kerr black hole. Roughly 
speaking, the measurement of one complex frequency (two observables) provides us with a determination of the mass 
and angular momentum of the black hole, and the measurement of the second frequency allows a consistency check 
with the general relativistic predictions [1, H, d, 0] ■ 
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The detectability of different QNMs depends on the relative QNM excitation, which in turn is determined by the 
cause of the perturbation (i.e. by the initial data). Given a detection, the resolvability of different modes depends 
on the nature of the multimode ringdown waveform. Two scenarios are possible U I ^ I' or m ^ m' (the modes 
have different angular dependence) the angular scalar product between the modes is zero to a good approximation 
Q, and the SNR can be expressed as a sum in quadrature of the single-mode SNRs. In Q we called these multimode 
ringdown waveforms quasi- orthonormal. If instead I = I' and m — m' , but n =^ n' (we look at different overtones with 
the same angular dependence) the angular scalar product is very close to unity, and mixed terms appear in the SNR 
calculation and in calculations of parameter estimation accuracy. We may call these waveforms quasi-parallel. 

A quantitative estimate of the relative QNM excitation for quasi-orthonormal waveforms will require numerical 
relativistic simulations and realistic initial data for the merger However, for quasi-parallel waveforms perturbation 
theory (combined with a good approximation to the initial data for ringdown) can provide useful information on the 
relative excitation of different overtones. In this paper, using perturbation theory, we develop tools to estimate the 
relative excitation of different overtones of Kerr black holes for general classes of initial data. We also try to gain 
some theoretical insight into the initial-data dependence of the results for realistic mergers, considering simple classes 
of initial data (such as localized spikes and gaussians) as model problems. 

Our work is the first application to gravitational perturbations of Kerr black holes of a theoretical framework 
developed over the years by different authors. Following the pioneering numerical analysis by Vishveshwara p^ . who 
studied the scattering of gaussian wave pulses on the Schwarzschild background, a number of studies investigated 
the analytical structure of the Green's function that propagates small perturbations in black hole geometries. Our 

but the Green's function in 



own analysis is based on the formalism developed by Leaver in the eightie s 111 
the Schwarzschild background has been studied by many authors |14| . [isl Il6l 



Il 7| . | 18| . More recently Glampedakis 

and Andersson extended the analysis to scalar perturbations of Kerr black holes [l9ll2C| . Here we carry out the first 
study of general-spin perturbations of Kerr black holes, including the physically most important case of gravitational 
perturbations. 

The main result emerging from [H [13, El, Q [S El, E [il E Hi is that a "universal" (initial-data independent) 
measure of the relative QNM excitation is encoded in the poles of the Green's function that propagates small pertur- 
bations of the geometry. These "universal" quasinormal excitation factors (QNEFs) are defined in Eq. (|3.8p below. 
They depend only on the Kerr geometry, not on the details of the perturbation. When combined with a knowledge 
of the initial data they can be used to compute the so-called quasinormal excitation coefficients defined in Eq. (j3.9p . 
which are a concrete measure of the QNM content of a waveform. 

The paper is organized as follows. In Sec.|lTl to develop some physical intuition, we consider a very simple physical 
system: a vibrating string with fixed ends. This part provides a useful pedagogical introduction to the QNM excitation 
problem, and may be skipped by readers familiar with the topic. We consider the Green's function solution of the 
vibrating string equation for generic initial data, identify the normal modes of the system as poles of the Green's 
function and show the importance of initial data to determine the excitation of the modes. In Sec. IIIII we show 
that many features of the vibrating string problem carry over to black hole perturbation theory, stressing the main 
differences between normal mode expansions and QNM expansions. We also anticipate some results on the convergence 
of the QNEFs, which are presented in more detail in Appendix El In Sec. lIVI we outline our calculation of the QNEFs, 
present the numerical results and compute the response of a Kerr black hole to localized and gaussian initial data. 
Technical details, as well as a discussion of the numerical accuracy of our calculations, are relegated to Appendix 
IB] Appendix [C] clarifies the relation between the Teukolsky and Sasaki-Nakamura (SN) formulations of the Kerr 
perturbation equations, and between the corresponding QNEFs. 

II. A PEDAGOGICAL EXAMPLE: A VIBRATING STRING 

Some key features of the black hole perturbation problem we address in this paper are exemplified by a very simple 
system: a vibrating string with fixed ends. To simplify the mathematics we pick units so that the velocity of the waves 
in the string c = 1, and consider a string of length tt. Then any disturbance of the string obeys the wave equation 

^ = ^ on 0<x<., (2.1) 
with u{t, 0) = u{t, tt) — 0. The general solution of this problem is easily verified to be 

oo 

u{t, x) = ^ [Cn COS nt -\- C'^ sinni) sin (nx) , (2.2) 

n=l 
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where we used an overbar to avoid confusion with the quasinormal excitation coefficients C„, as defined in Eq. (j3.9p 
below. In Fourier language we say the general solution is a superposition of normal modes with sinusoidal dependence 
on X and t, labeled by an integer n. Each mode has frequency ui = n. 

For the general solution to be useful we must determine the constants C„ and C^, that is, we must determine the 
contribution of each individual mode. This contribution can easily be computed once we are given initial data, namely 
the initial configuration u(0,x) = uq{x) and velocity profile dtu(0,x) = vo{x) of the string. Indeed, consider (12. 2p and 
its first derivative, both evaluated at t = 0. Multiplying both sides by sinnx and integrating on (0 , tt) we get 

Cn — — uo(x) sinnxdx, C' — — / vo(x) sinnxdx, (2-3) 

71" Jo Jo 

which completely specifies the solution. 

In more general situations it is not possible to find closed-form elementary solutions satisfying some given boundary 
conditions (in our vibrating string example, the fixed ends condition). However, an elegant formal solution can be 
obtained using Green's functions. Let us consider a shght generalization of Eq. (|2.ip : 

(2-4) 

where we introduced a potential V{x) and a source S representing, say, external forces acting on the system (for a 
free vibrating string V{x) =5 = 0). We define the Laplace transform of u(t,x) as^ 



Cu{t,x) = u{uj,x) ^ u{t,x)e"^'dt, (2.5) 



In terms of the Laplace transform, the original field can be written as 



1 



oo+'ic 



u{t,x) = — u{uj,x)e-'"^*duj. (2.6) 

Using the elementary property C ^"g*'^-* = —iujCu — e*"*°u(a;, to), the Laplace transformation of (|2.4p leads to 

^+[u:^-V{x)\u^I{lo,x), (2.7) 

where 

du{t, x) 



iuju(t, x) 



dt 



I{lo , x) = e*"*" 

t=t 

This equation is formally solved with the use of a Green's function G{x,x') such that 



S. (2.8) 



dx'^ 

In terms of the Green's function the solution is simply given by 



" + [cj2 _ V{x)\ G{x , x') = 5{x - x') . (2.9) 



I{uj,x')G{x,x')dx'. (2.10) 



Suppose we know the Green's function. Then the previous equation shows that, given I(lu ,x') (which means, in the 
absence of external forces, given initial data) we can determine, at least in principle, the solution. 

There is a general prescription to construct the Green's function (2T| . Find two linearly independent solutions of 
the homogeneous equation, say ui(uj ,x) and U2i^,x), each satisfying one of the required boundary conditions: for 
the vibrating string these solutions would be such that -1*1(0; , 0) = 0, U2{'-^ 7 ^r) = 0. The Green's function is then 



1 1 {ii{x)u2{x') if a; < x' , 
G{x,x') = — { ^ ^ ^ (2.11) 

{ll{x'){l2{x) a x' < X , 



w 



^ The usual Laplace variable s = —iuj. We prefer to use ui for notational consistency with previous work by Leaver [l3j and Andersson 
[iTll . Our transform is well defined as long as Im(a;) > c. 
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where we defined the Wronskian between the two solutions 



dill . . du2 
W = -t:—U2 - ui—— , 
ox ox 



(2.12) 



which for equations of the type (|2.7p is a constant. For the vibrating string the homogeneous solutions are elementary 
functions: ui = sinux, 112 = sina;(a; — tt), and the Green's function 



\ujxs\nuj[x' — tt) \^ X ^ x' 



G{x,x') = 



UJ sm UJ7T 



(2.13) 



UJ sm UJTT 



if x' < 



Notice that the Wronskian W — —wsmajn is zero at w = n with n integer, that is, at the normal frequencies of the 
system. Near the poles uj = n we have W ~ — [n7rcos(n7r)] {uj — n). Corresponding to zeros of the Wronskian (which 
are also poles of the Green's function: see below) the two solutions iii and U2 are no longer independent: they satisfy 
both boundary conditions simultaneously. In fact, setting lu — n in (|2.13p we can see that ui and 112 coincide, and 
correspond to the normal modes of the system. Using (|2.6p . (|2.8p and (|2.10p and setting for simplicity to = S = we 
get 



u{t,x)^^ I dx'duj [icuuoix') - vo{x')] G{x , x')e-"^* 

ZTT 



(2.14) 



The w— integral can be performed by closing the contour of integration. We choose the contour depicted in Fig. [TJ 



— X — X — X — X — X — n — X — X — X — )f- 



^ — )^ 

Rc(^) ;' 



FIG. 1: Integration contour for the vibrating string problem. Crosses mark zeros of the Wronskian W , corresponding to the 
normal frequencies of the system. 

If ui{x) , U2{x) are analytic and have no essential singularities inside the contour^, the poles of the Green's function 
are all due to zeros of the Wronskian. In this way we get for the integral in (|2.14p . 

u{t,x) = — [ dx'duj[iujuo{x') ~ vo{x')]G{x ,x')e-"^* 
27r J 

= —i'S^ — i / dx' sinnx' [inuo{x') — vo(x')]\ e~^"'* sinnx 
^ rnr IJ J 

duj / dx' [iujUQ{x') ~ vo{x')]G{x ,x')e~^'^* , (2-15) 



1 

'2^ 



HC 



^ Both conditions are met in the case of a vibrating string, but both are violated when we deal with black hole QNMs. For black holes, 
the Green's function essential singularity at the origin gives rise to tails, and the integral over the half-circle at infinity is responsible 
for the early-time response of the black hole. See discussion after Eq. 1 13. 6I I. 
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where "HC" in the second term means that the integration should be performed along the half-circle. Not surprisingly, 
taking the real part of the right hand side we recover the result (|2.2p with expansion coefficients given by (|2.3p . 

The bottom line of this discussion is that the general solution depends crucially on two elements: (i) the residues 
of the Green's fimction evaluated at the poles (that is, at the normal frequencies); (ii) the function /(a;,x), which (in 
the absence of initial forces) is nothing but the initial data. By inspection, the net result for the field can be expressed 
as a sum of the form 

where = n is a normal frequency, J„(t ,x) — J dx' I{ujn , x') Un{ujn , x') and u„(cl;„ , x) = sinna; is a normal mode 
wavefunction (that is, any of the homogeneous solutions evaluated at the normal frequency aj„ n). We will see 
below that a similar result holds for the ringdown of Kerr black holes. 



A. Effect of initial data 



In the previous section we pointed out that initial data play a crucial role to determine the excitation of the normal 
modes of a system. For illustration, below we consider three simple examples that will be useful in the following to 
understand, by analogy, the initial-data dependence of the excitation of a Kerr black hole. 

Suppose first that we have an (initially stationary) plucked string: the string's initial profile is a triangle of height 
h with a vertex at x^ ^ i.e. vq{x) = and 




\iQ<x<x^ , 

(2.17) 



— it X < a; < TT 



Stationarity of the initial data implies that C'j = 0, so that 



LKt^x) — > C„cosntsmnx, = 5- 5^ (2.18) 

^-^ ^•'(vr — x°) 



n=l 



Notice that modes having a node at x^' , where the string is plucked, are not excited ((?„ — 0). Notice also that the 
"excitation factors" Cn decrease as 1/n^. 

As a second example take stationary, localized "(5-function" initial data of the form 

m{x)^S{x~x'^), V(){x)=0. (2.19) 

The excitation factors are trivially computed: 

Cn = - sin {nx'^) . (2.20) 

TT 

It is apparent that all modes are excited to a comparable amplitude except for modes with a node at the plucking 
point, which are not excited at all. 

Our third and last example are stationary, gaussian initial data: 

uo(^) - , vo{x)=0. (2.21) 

For large b the gaussian is strongly peaked at a; = x"^, in which case the contribution to the integral outside of (0,7r) 
can be ignored and we have 

Cn = - r uo{x) smnxdx ~ - T e-H— sinnxdx = ^ sin(nxg) ^.„V(4fc) ^2.22) 

Therefore a mode with given n is maximally excited when the width of the gaussian satisfies the condition b = r? 12. 
A similar result will be seen to hold for gaussians exciting Kerr black holes. 

The basic lesson we learn from these examples is that the excitation of a system is very sensitive to the initial data. 
More specifically, whether a given mode is excited or not depends strongly on the point where we excite the system 
("pluck the string"). 
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III. OSCILLATING BLACK HOLES 



In the SN formalism [22, [23, [2J, [2^ , perturbations of a Kerr black hole induced by a spin-s field are described by 
a single function X(''^(t, r) whose Laplace transform satisfies 

^^^^^ + VsnX^-^^ , r) ^ /(c , r) , (3. 1) 

where the effective potential Vsn depends both on the radial coordinate and on the frequency. The function I(uj ,r) 
is a linear combination of the SN function X'-^^{to , r) and its time derivative X^^\tf) , r) at time to (see below for the 
explicit expression for scalar perturbations). The tortoise coordinate is defined by the condition 



dr 



(3.2) 



ranging from —oo (the location of the event horizon) to +oo (spatial infinity). We use Boyer-Lindquist coordinates 
and follow Leaver's choice of units, setting G = c = 2M = 1. In Leaver's units the angular momentum per unit mass 
a is such that < a < M = 1/2, and the horizon function A = — r + a^. Sometimes we will present our results 
in terms of the more familiar dimensionless angular momentum j — 2a, such that < j < 1- The class of problems 
that fit in this description include any massless field in the Kerr geometry, including gravitational, electromagnetic 
and scalar fields [2^ . 

The main difference with the vibrating string example is that our system is not conservative: waves can escape to 
infinity. For this reason an expansion in normal modes is not possible (see [l3, [iH [13] for extensive discussions of this 
point). Wave pro pagation is also complicated by backscattering off the background curvature, which is responsible 
for tail effects (28|. Despite these complications it can be shown that the poles of the Green's function (now located 
at complex frequencies corresponding to the QNMs) still play an important role in the evolution. 

The QNM contribution can be isolated from other features of the signal, such as the late-time tail, using the Green's 
function technique (isl . [l7j . First one defines a solution of the homogeneous equation having the correct behavior at 
the horizon (only in-going waves), 

lim Xj.'^ - ^~^{^-mn)r, ^ ^3 3) 

lim X(^) ~ AMe-'""^' + Aout(o7)e'"'^* , (3.4) 

where r± — [1 ± (1 — 4a^)^/'^]/2, and a second solution x!^\ behaving as e'"''* for large values of r. Since the 
Wronskian W — 2iujAin we can express the general solution as [l3] 

To proceed we make the astrophysically reasonable assumption that the observer is located far away from the black 
hole. If the initial data have compact support and this support is entirely located closer to the black hole with respect 
to the observer (this is basically a "no- incoming radiation from infinity" condition), a good approximation will be 

X^-'\^ ' ^ 2kJA~ ^("^ ' '^^rl^'^K • (3.6) 

As explained in (isl . [itI [29j , when we invert this expression to get the solution in the time domain we get (as in the 
vibrating string case) a contribution from the poles of the Green's function. This contribution can again be isolated 
by closing the path of integration, as was done in Eq. (j2.15|) . An important difference that distinguishes black hole 
spacetimes is that there is now an essential singularity at = and a branch cut extending from the singularity 
to — ioo (29l. [30I [Slj. To prevent the essential singularity from lying inside the integration contour we must modify 
slightly the contour of Fig. [T] We place a branch cut along the negative imaginary-w axis and split the half circle at 
00 into two quarter circles. The new contour is shown in Fig. [2l The poles in the complex frequency plane 
are the zeros of ^in: they correspond to perturbations satisfying both in-going wave conditions at the horizon and 
out-going wave conditions at infinity, that is (by definition) to QNMs. QNM frequencies have negative imaginary 
part. Since we assume a Fourier dependence of the form X'-'^'>{t, r) ~ ^-^^t this means that there are no exponentially 
growing modes. There is an infinity of QNMs, which are usually sorted by the magnitude of their imaginary part and 
labelled by an integer n. 
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FIG. 2: Integration contour to invert Eq. (|3.6|l . The shaded area is the branch cut and crosses mark zeros of the Wronskian W 
(the QNM frequencies). 



Extracting the QNM contribution^ to the radiated wave we get 



X'^'\t,r) = -Re 



An 



-Re 



E 



(3.7a) 



(3.7b) 



where the sum is over all poles in the complex plane and the quasinormal excitation factors (QNEFs) _B„ are defined 
as 



Bn 



Ao 



A-out f dA{, 



2Luan 2lu \ duj 



(3.8) 



Here a„ is a commonly used notation for the derivative of at the QNM frequency. The quasinormal excitation 
coefficients Cn can be computed as 



Cri — Bn 



(3.9) 



whenever the integral on the right hand side, which must be evaluated at the QNM frequency co — w„, is convergent. 
In general the QNM frequencies w„, the B^s and the C„'s (as well as the wavefunction) depend on {l,m) and the 
spin of the perturbing field s, but to simplify the notation we will omit this dependence whenever there is no risk of 
confusion. 

By definition of the QNM frequencies, Xr'^^ ~ Aoute'"'^* as ^ cxd at = w„: in this sense the above integral 
is "normalized". For source terms I{uj ,r) that are zero outside some finite range of r, or have sufficiently rapid 
exponential decay as |r*| oo, the integral is also convergent. For other classes of initial data the integral as 
evaluated on the real line is, in general, divergent. This is a major difference with respect to ordinary normal mode 
expansions. The normal modes of closed mechanical systems are Sturm-Liouville eigenfunctions of the wave equation, 
and their excitation coefficients are weighted integrals of the source term over the mode. Since the Sturm-Liouville 
eigenfunctions are always bounded, the integrals always converge. For QNM expansions a meaningful definition of 



^ There are also other contributions to the signal. The essential singularity at the origin (a; = 0) gives rise to the tail of the time evolved 
wavefunction; the integral over the quarter circles at infinite frequency produces the early time response of the black hole. Since our 
main focus is the QNM contribution we discard these two terms in the integral. For more details we refer the reader to [I3l |29| . l3l | . 
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the integral for general source terms requires more care, and the introduction of an analytical continuation procedure 

MM- 

Equation (|3.7p is one of the main results we will use throughout the rest of the paper. Once we specify initial 
data, the QNEFs allow the determination of the QNM content of a signal. QNEFs have long been known for scalar, 
electromagnetic and gravitational perturbations of Schwarzschild black holes [ij, [13, ITtI. More recently there have 
been attempts to extend those calculations to scalar perturbations of Kerr black holes [id [20l | . Electromagnetic and 
(most importantly for gravitational-wave phenomenology) gravitational perturbations of Kerr black holes have not 
been dealt with so far. One purpose of the present work is to fill this gap. Before turning to an explicit calculation 
of the BnS we will provide the explicit form of the function /(w , r) for s = 0. We will also address two important 
conceptual issues: the time-shift problem and the convergence of the QNM expansion. 



A. Initial data for scalar perturbations 

For a scalar field fl^ the function I{uj , r) is given by 

2amr + uJO^Aji 



/(a;,r) 



(r2 



1 



where 7;,„ is a number depending on the multipolc of the field under study: 



7Jn 



27r 



Si^{auj, d)Simiauj, 9) sh 



(3.10) 



(3.11) 



Here the Sim's are scalar (s = 0) spin- weighted spheroidal harmonics Q. We require them to satisfy the normalization 
condition 27r Sf^ (aui, 6)Sim. (auj, 9) sin 9d9 = 1 (notice that this normalization differs by a factor 27r from that 
adopted in [l9f). Even though these functions depend on aw, in Q we showed that in some cases (for example, when 
we are dealing with slowly damped QNMs) this dependence can safely be neglected. When aoj — scalar spheroidal 
harmonics reduce to the usual spherical harmonics, and the integral can be computed analytically with the result 



"fin 



2{P 



(3.12) 



For the calculations in this paper, the analytic formula is extremely accurate even when atu ^ 0. For example, for 
a near-extremal black hole with j — 0.98 a numerical evaluation of the integral at the fundamental scalar QNM 
frequency yields 722 = 0.854139, in excellent agreement with the analytic prediction 722 = 6/7 ~ 0.857143. Results 
for other slowly-damped modes are similar. 

In the following of the paper we will focus, for simplicity, on static initial data. For large r Eq. (|3.10p reduces to 



/(w,r) ~«a;e*'^*''X(°)(io,r), 

whereas for r very near the horizon we get 

/(w , r) = ie*"*" (uj - 2m£t) X^^^ {to , r) 
where = a/r+ is the angular velocity of the horizon. 



(3.13) 
(3.14) 



B. The time-shift problem 

Taken at face value Eq. (|3.7p is troublesome, because the response diverges (exponentially) for early times. The 
well-known fact that ringdown waveforms only make sense for certain values of t is usually referred to as the "time-shift 
problem" MM, 

[1^. For illustration, consider Eq. (|3.7p for a Schwarzschild black hole with static initial data: 



-Rc 



n=0 



-to) 



Ao 



dr' 



(3.15) 



The wave field X('*)(t,r) at some observation point r and time t depends exponentially on the (arbitrary) choice 
of to J which does not make much physical sense. Suppose in addition that the initial data consist of a narrow pulse. 
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say for simplicity a J-function: X'^''\to,r) = 6(r — r^). Far away from the black hole Xr^ ~ ^□ute*""''*, and the 
response of the black hole increases exponentially with . This is against physical intuition: according to our Fourier 
convention stable oscillations have Im(a;) < 0, so a small bump in initial data far from the black hole seems to excite 
much more ringing than a huge bump in the strong-field region around the horizon. 

The way out of this problem [14 . TSj is to realize that ringdown only starts after the initial data (that for simplicity 
are usually assumed to have compact support and be located in the far zone, ^ 1) reach the potential barrier, 
where ringdown originates and the perturbation travels back to the observer. Suppose the observer is sitting at 
some large and an impulse localized at r* = rf falls into the black hole at time to- The Regge- Wheeler "potential 
barrier" has a location depending on the mode under consideration, and the very notion of potential barrier becomes 
fuzzy for Kerr black holes. A reasonable estimate is to assume that the "barrier" is located close to ~ 0, and 
to avoid complications we will take this estimate as our fiducial value. In geometrical units, based on the above 
discussion, the starting time for the ringdown signal will be approximately 

tstart = |rf|+n +to. (3.16) 

For the QNM expansion to make sense, the observation time t must always be chosen so that t > tgtart- This choice 
partially gets rid of the time-shift problem, in the sense that to and rf do not appear explicitly in the waveform. 

A more rigorous approach should also take into account the contribution of the prompt response (the analogous of 
light cone propagation in flat space) and the matching between these two phases. Ref. [l8| explores some aspects of this 
problem, computing the large-frequency contribution to the Green's function (large frequency is roughly equivalent 
to early times). The final result does not significantly alter our "simple-minded" scenario. We should also mention 
that our prescription to solve the time-shift problem works reasonably well only for sharply localized initial data, 
like narrow gaussians. If the initial data have a significant spread in the radial direction the starting time tgtart is 
ill-defined and we must resort to a more detailed analysis, based perhaps on time-dependent excitation coefficients 
p^ . This problem deserves further investigation. 



C. Convergence of the QNM expansion 

The convergence of the QNM series has been studied in some special important cases by Leaver . In particular 
he pointed out that the sum can only be expected to converge at times t such that ringdown dominates (see the 
above discussion of the time-shift problem), and that the ultimate convergence of the series depends on the large-n 
asymptotics of the QNEFs. The convergence problem was revisited numerically by Andersson computing the 
Bn's up to n ^ 200 for scalar perturbations. His calculation shows that their magnitude decreases monotonically for 
large n. From (|3.15p it follows that, if the initial data are localized at r* = rf ^ 1, the ratio of two consecutive 
terms behaves like (Bn+i/ Bn)e~^'-^"+^~'^"'>'-*~^*~^*\ Asymptotically, LUn+i ~ — i/2 for large n. In Appendix [A1 
we determine analytically, for the first time, the asymptotic behavior of the QNEFs for a number of static black hole 
spacetimes. Our most interesting result is that the large-n behavior of scalar and gravitational perturbations of a 
Schwarzschild black hole is 

= 4-2) = - \ , (n oo) , (3.17) 

3 [± log 3 — (2n + l)7rzj 

in good agreement with numerical calculations (cf. Figure 2 in Ref. [iB|). This implies that \Bn+i/Bn\ ^ 1 in 
the same limit, so the ratio of consecutive terms is of order exp[— (t — r^, — )/2] and the sum will converge for 
(t — r» — rf ) ^ 0. This is precisely the kind of convergence we should expect from a QNM expansion. 

To prove more rigorously the convergence of the expansion (j3.7p we still need two ingredients. The first is the 
convergence of the integral in (j3.9p . For generic initial data this integral diverges on the real line. The divergence can 
be cured by analytic continuation, evaluating the integral as a contour integral (for details see [13, [l3|)- Secondly, the 
fact that the QNEFs _B„ ~ 1/n for large n does not imply convergence of the corresponding quasinormal excitation 
coefficients C„ . A proof of convergence would require the calculation of the integral over initial data and a knowledge 
(or an estimate of) its n-dependence for large n. 

Having discussed these basic properties of QNM expansions, in the following of the paper we turn to an explicit 
calculation of the j-dependence of the QNEFs for different QNMs and for different spins of the perturbing field. 
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IV. QUASINORMAL EXCITATION FACTORS IN THE KERR BACKGROUND 

A. Formalism 

Perturbations in the Kerr geometry can be studied using both the Teukolsky [35| and SN formahsms. There 
are two main advantages in using the SN formahsm: (i) the potential and source in the SN wave equation are short 
ranged and (ii) in the hmit a — > the SN wave equation reduces to the Regge- Wheeler equation, describing axial 
perturbations of a Schwarzschild black hole. 

The QNEFs Bn depend, by their own definition ()3.8|1 . on the amplitude of ingoing and outgoing waves at infinity. 
Even in the Schwarzschild limit, despite the isospectrality between axial and polar perturbations, the Regge- Wheeler 
and Zerilli QNEFs differ by a proportionality constant. Denoting quantities in the Zerilli equation by a plus (meaning 
even-parity), the two are related in a simple way IS]: 

^ {l-l)l{l + l){l + 2) + 6iu; (/-l)^a + l)(/ + 2) + 6»^ 

{I - l)l{l + + 2) - 6iiu """^ ' " {l-l)l{l + l){l + 2)-6iuj ^ ' 

Similarly, the SN QNEFs for spin-s perturbations S'*^ will differ from the Teukolsky QNEFs (that we will denote 
by a subscript T, B^^) by a proportionality constant. In Appendix [Cl we derive the following relations: 

4"'=i?(°), 4-'^=-:^ ^ ^S^-'^ 4-^^^ ,,, , (4.2) 

zamuj — Aim ~ a uj'^ A(A + 2) — oilu — 12auj(aLU — mj 

with A = Aim + (aw)'^ — 2amu!. 

Our calculations in this paper rely heavily on the formalism developed by Leaver [ill, fisl ] , which refers to the 
Teukolsky formulation of the Kerr perturbation equations. In this Section we only outline the basic steps of the 
calculation, relegating details to Appendix [Bl Expanding a spin-s field ijj{t ,r ,9 ,(/)) as 

t^(t,r,e,(f>) = — I e-^'^'Yl ^ e""^^imWi?imWdu;, (4.3) 

l — \s\ m— — l 

we get ordinary differential equations for Sim and Rim [33 |: 

du [(1 - u^)duSim] + f a^u?u^ - 2aujsu + s + - \ Sim = , (4.4) 

Ad^Rim + {s + l)(2r - l)drRlm + VRim = , (4.5) 

where u = cos 9 and 

V = 2iswr - a^ijj'^ - Aim + [(^"^ + a^ fuj'^ - 2amujr + a^m^ + is {am{2r - 1) - w(r^ - a^))] . (4.6) 

The field spin s = , — 1 , —2 for scalar, electromagnetic and gravitational perturbations, respectively. The separation 
constant Aim reduces to l{l -I- 1) — s{s + 1) in the Schwarzschild limit. The general series solution of the angular 
equation with appropriate boundary conditions is given in Appendix IB II 

Let us consider the radial equation (|4.5p . For brevity, in the rest of the paper we drop the dependence of the 
Teukolsky function Rim on the angular indices {l,m). Following Leaver (see Appendix IB] for details) we introduce 
three different solutions of the homogeneous equation (|4.5p : ,Roc+ ,-Roo_, with asymptotic behavior [l^ 

lim Rr. - (r+ - r^)-^-'+''^+"'+e"^''+{r - r+)-'-"'+ , (4.7a) 

lim Rr^ - Al^{Lj)r~'^-''^e-"^'' + Al^^{cj)r-^-^'+"^ e''^'' , (4.7b) 



lim Roo. ~ r-i-2^+*"e*"'' , (4.8a) 
lim Roo_ r.-i-''^e~''^'' . (4.8b) 



In these relations (T-(- = {tor-^- — am) /b and 6 — vl — Aa^. 
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By definition, QNM frequencies w„ are sucfi that Af^{u!n) = 0. We compute tfiem following the procedure in [lll.l33j 
and verifying that Af^ = to some required numerical accuracy (typically better than one part in 10^). From (|4.7bp 
and (|4.8bp it follows that 

Rr+ = Al^Roo_ + Al^^Ro^^ , (4.9a) 
K+ = ^hi-R'oo_ + ^out^!x.+ • (4.9b) 

From these relations it is clear that computing the QNEFs requires accurate representations of the solutions both 
close to the horizon and at infinity. Close to the horizon, an accurate representation is provided by a series solution 
first found by Jaffe. At infinity we must resort to a different expansion in terms of Coulomb wavefunctions. Properties 
of these representations which are important for the calculation of the three solutions (-Rr+, -^00+ and i?oo_) are given 
in Appendix IB 21 and IB 3[ respectively. We performed a number of consistency checks on the solutions, some of which 
are described in Appendix IB 41 More details can be found in the original papers by Leaver [ll], [l^, [l3| • Notice that 
the outgoing- wave amplitude Aout is unambiguously defined once we fix the normalization according to Eq. (|4.7p . The 
corresponding normalization for the SN wavefunctions is derived in Appendix [Cl see in particular Eqs. \C7\ . (jC13p 
and (lUIel) . 

To compute the QNEFs (14. 2p as functions of the rotation rate j — 2a, we first obtain the amplitudes AT^ and 
A^^^ by matching the three solutions and their derivatives at some point (usually r = 5 — 7). Then we compute 
the derivative a„ introduced in (|3.8p . Close to a QNM frequency we can perform a Taylor expansion of the ingoing 
amplitude: ~ a„(5aj, with Soj a small complex number. For each value of s and a we evaluate Af^ at a discrete 
set of points co = uin + kSu (typically we pick Su) = 10~* and the integer k = —2, . . . , 2). Then we obtain a„ from a 
linear fit of Af^ as a function of 6u!, and verify that the resulting derivatives satisfy the Caucliy-Riemann conditions. 

In the rest of this paper we present our results for the QNEFs using both the Teukolsky and SN formalisms, and 
discuss some of their implications for gravitational-wave phenomenology. 

B. Numerical results 

The QNEFs of Schwarzschild black holes were first introduced and computed by Leaver in his seminal analysis 
of the radiative Green's function. Table I of [l^ lists the QNEFs for the first four overtones of electromagnetic 
perturbations with 1 = 1, the first seven overtones of gravitational perturbations with / = 2 and the first four overtones 
of gravitational perturbations with I = 3 and I = 4. Building on Leaver's analysis. Sun and Price [l3| studied the 
initial data dependence of the Schwarzschild quasinormal excitation coefficients. Andersson computed the QNEFs 
for scalar and gravitational perturbations using the approximate phase-integral method, finding good quantitative 
agreement with Leaver for s = —2 and pointing out a few sign errors in Leaver's results. Later, using the same 
technique he was able to compute the (scalar) QNEFs up to overtone numbers n ^ 200 {l^- Glampedakis and 
Andersson |20| tried to extend the methods of Ref. |l3| to scalar perturbations of Kerr black holes. The QNEFs listed 
in their Table 3 do not reduce to the correct Schwarzschild limit [s^], and the "effective amplitude" results in Figure 
5 of their paper are only qualitatively correct. 

We computed QNEFs using both Mathematica and a Fortran code. In the limit a ^ the Fortran code reproduces 
Table I of pJi], once we correct for Leaver's sign mistakes, to all (five) significant digits. Given the complexity of the 
algorithm described in Appendix[B]it is hard to quantify our numerical error for a 7^ 0, but results should be accurate 
to the same level. The numerical accuracy depends on many factors, including the number of terms included in the 
Coulomb wavefunction expansion (|B39p . the matching radius chosen to solve the linear system (|4.9I) and the value of 
6u! used to compute a„ by a linear fit. In general, the number of digits to which ylin = at the QNM frequency and 
the level to which the Cauchy-Riemann conditions are satisfied in the calculation of a„ are good accuracy indicators. 
We decided to list the QNEFs with a five-digit accuracy, but in some cases (especially for corotating modes with 
I = m and large rotation rates) the number of significant digits may be smaller. 

1. Scalar perturbations 

In Tables HI and [III we provide the QNM frequencies and damping times for scalar perturbations with I = 2, all 
values of m and selected values of the angular momentum j = 2a. In Table [1111 we list the corresponding QNEFs, 
and in Table II VI the outgoing-wave amplitudes Aout- In the Schwarzschild limit, our QNEFs can be compared with 
Andersson's |l7l. [l8|. His result for the fundamental scalar mode with I = 2 agrees with ours to a five-digit accuracy 
[35| . This agreement is quite impressive, given the approximate nature of the phase- integral method. For / = his 

result for the fundamental mode (-Bq"'' — 0.197 — 0.046i) is within ~ 7% (~ 23%) of the real (imaginary) part of our 
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result: Bq ' ~ 0.212 — 0.059i. The agreement is still very good, considering that the phase-integral method is based 
on a WKB-type approximation and should only be accurate in the "eikonal" (large I) limit. 




FIG. 3: Real part (left) and imaginary part (right) of the scalar QNEFs for the fundamental mode with I = 2. Solid lines 
correspond to |m| = 2, dashed lines to |m| = 1, the dotted line to m = 0. The lines displaying oscillations for large j (black 
and green in the colour version) have m > 0. 



Some physical insight can be obtained by plotting the real and imaginary parts of the scalar QNEFs for the 
fundamental mode with ^ = 2 as functions of j (Fig. [3]). The excitation of modes with m < is a slowly- varying 
function of the rotation rate. On the contrary, as j ^ 1 the real and imaginary parts of -B^"^ for corotating modes 
become rapidly oscillating functions of j. Our calculation is in remarkable agreement with a classic result by Ferrari 
and Mashhoon |36|] : when I = m the QNM excitation tends to zero as j ^ 1 , and this is an indication (if not a proof, 
given the incompleteness of QNMs) that extremal Kerr black holes are not marginally unstable. 



T ' 1 ' 1 ' 1 ' r 




Re(5) Re(5) 



FIG. 4: Path of the scalar QNEFs for the fundamental mode (left) and first overtone (right) with I = 2 as, j varies in the range 
< j < 0.996. Linestyles are the same as in Fig. [S] 

More features can be seen in Fig. IH where we plot the path followed by the real and imaginary parts of Bn\ 
thought of as parametric functions of j, for < j < 0.996. The oscillations of the real and imaginary parts of B*^"' 
for the corotating mode with I = m = 2 produce a spiral in the complex plane. For small j, the QNEFs of corotating 
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and counterrotating modes with the same |m| move in opposite directions, approaching the Schwarzschild limit with 
the same tangent. Corotating and counterrotating QNM frequencies tend to their Schwarzschild limit in a similar 
way [33I. [37}. Another noteworthy feature of Fig. [4] is the scale of the real and imaginary axes: for corotating modes 

with I = m and large rotation, |-BJ''^| is roughly five times larger than I-Bq"^!- other words, for perturbations with 
I = m the contribution of the overtones should be more significant for rapidly rotating black holes. In the following 
we will see that many of these considerations are still valid when we consider perturbations of spin s ^ 0. 

2. Electromagnetic perturbations 

It is widely believed that astrophysical black holes should possess very little charge (if any). For this reason 
electromagnetic perturbations of a Kerr black hole (s = —1) arc not considered of great astrophysical relevance. 
Nonetheless, electromagnetic perturbations are of more than academic interest. For example, they could find useful 
applications in models of elementary particles based on the Kerr-Newman metric (see (38j] and references therein). 




FIG. 5: Path of the electromagnetic QNEFs for the fundamental mode (left) and first overtone (right) with I = 1 as j varies in 
the range < j < 0.996. The top panels refer to the SN formalism, the bottom panels to the Teukolsky formalism. Linestyles 
are the same as in Fig. [S] 

Results for the lowest radiative multipole I = 1 are presented in Fig. [5l Compared with the scalar case, the 
main new feature is that QNEFs differ (both analytically and numerically) depending on whether we use the SN or 
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Teukolsky formalisms: see Eq. (|4.2p and Appendix [Cl Of course, both formalisms are equally legitimate, and the 
use of one or the other depends on the physical problem at hand. By construction the SN QNEFs reduce to the 
Regge- Wheeler QNEFs in the Schwarzschild limit, and in this sense they have a more direct physical interpretation. 
As j ^ the (Bardeen-Press-)Teukolsky quantities can be transformed to the corresponding metric quantities in the 
Zerilli (Regge- Wheeler) formalism using the differential transformations derived by Chandrasekhar [39|. 

Apart from this distinction, the qualitative features of Fig. [5]are similar to the scalar case of Fig. |4l For small j, the 
QNEFs of corotating and counterrotating modes with the same |m| approach the Schwarzschild limit with the same 

tangent. For corotating modes with I = m = 1 and large rotation rates, \b[ ^''| is significantly larger than \Bq 
This is another indication that, independently of the the value of s, the high-overtone contribution should be more 
significant for I = m modes and rapidly rotating black holes. 



3. Gravitational perturbations 

Gravitational QNEFs in the Teukolsky and SN formalisms are listed in Tables IVllVIl for the first two overtones with 
I = 2, and selected values of the black hole's angular momentum. The corresponding QNM frequencies and angular 
separation constants can be found in Tables II and V of [4], respectively (there we list the imaginary parts with the 
opposite sign). As j our numbers agree to all digits with Table I of [l^, except for the relative sign of a few 
modes with I = 3 and I = 4. These minor sign errors in Leaver's paper were first pointed out in [l7j . 

In Fig. [6] we plot the real and imaginary parts of the gravitational QNEFs for the first two overtones with / = 2 in 
the SN formalism as functions of j. The qualitative behavior is remarkably similar to the scalar case (cf. Fig. [3]). For 
fast rotation, the maximum excitation of the first gravitational overtone (as compared to the maximum excitation of 
the fundamental mode) is larger by a factor ^ 8. 

In Fig. [7] we plot the path followed by the real and imaginary parts of i?i "^^ , thought of as parametric functions 
of j, in the range < j < 0.996. There's a remarkable similarity between the plot for n = in the SN formalism 
and the scalar case of Fig. ID Even though it's hard to accurately compute the QNEFs when j is very close to one, 
our numerics show clear evidence that in both cases the excitation of I — m modes is zero in this limit: the center of 
the spiral is located at the origin of the complex plane. Our calculation provides further evidence that extremal Kerr 
black holes are stable, and that QNMs of a fast rotating Kerr black hole are very hard to excite [l9l. [20l. [s^ . 

In Fig. [8] we plot the modulus of gravitational and scalar QNEFs for the fundamental mode and first overtone 
with I = m = 2. This plot is useful to visualize the effect of rotation on the excitation of different overtones, 
at least for radially localized initial data (for gaussian initial data the response is exponentially modulated by the 
QNM frequencies, see below). For gravitational perturbations the first overtone is maximally excited (with a relative 
amplification of a factor ~ 8 with respect to the fundamental mode) for j ~ 0.95. For scalar perturbations the effect 
is similar, but the relative amplification of the first overtone is smaller (roughly a factor 6) and the maximum occurs 
at j ~ 0.85. 



C. A simple application: localized and gaussian initial data 



As an application of our calculation of the QNEFs, we study the response of a Kerr black hole to initial data in 
simple situations amenable to an analytic treatment. The relevance of these model problems to realistic perturbations 
of astrophysical black holes is questionable, but our analysis could provide some insight into the black hole's response 
to more generic initial data. For simplicity in the following discussion we focus on scalar perturbations, so that 
Eq. (|3.10p applies. We start by considering static, localized initial data of the form 



X(°)(io,r) =(5(r,-rf), X '^"\t„ , r) ^ . 
Combining Eq. (|3.7p with the definition (|3.16p yields 



-Re 



^out 



2amr + w„a^7imA 
(r^ -f 0^)2 



For initial data in the far zone this is further simplified to 



-Re 



-iw„(t-t5tart) 



+00 . 



(4.10) 



(4.11) 



(4.12) 
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FIG. 6: Real part (left) and imaginary part (right) of the gravitational QNEFs with I — 2 and different values of the overtone 
index n, indicated in the inset. All plots refer to the SN formalism. Linestyles are the same as in Fig. |3] 



For initial data localized near the horizon, using (|3.3p we get 

rf^-c». (4.13) 

Rotational effects appear (in the form of the usual Lense-Thirring frame dragging) only when the 5-like disturbance is 
located near the horizon. Eq. (|3.17|) shows that in the Schwarzschild case -B„ ~ 1/n for large n. Since we expect the 
same leading-order behavior to hold for Kerr black holes, and we know that Kerr QNM frequencies scale as w„ ~ —in 
for large n |33'|, a ^-like source excites all modes to comparable amplitude. This is reminiscent of the analogous result 
(I2.20p for the vibrating string. 

The vibrating string analogy suggests that we may be able to associate nodes to the QNM eigenfunctions. In other 
words, there could be special locations of the initial data rf for which the black hole's response is zero (or at least 
very small) . We expect it should be easier to find these "nodes" , if they exist at all, for QNM eigenfunctions which 
are as close as possible to ordinary normal mode eigenfunctions. This happens when Re(a;„) S> Im(cj„), that is, for 
corotating modes with / = m and near-extremal (j ~ 1) black holes [s^l- For this reason below we focus on scalar 
perturbations of Kerr black holes with j = 0.98 and consider only modes with I = m = 1 ot I ^ rii ^ 2. 

In Fig. [5] we plot the response of a Kerr black hole with j — 0.98 (so that the horizon radius r+ 0.5994) to 
initial data localized at different values of , using Eq. (|4.1ip . We compute the QNM frequency tj„ and the angular 



X(°)(i,r) ~ -Re 



E 

n=0 
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FIG. 7: Path of the gravitational QNEFs for the fundamental mode (left) and first overtone (right) with I — 2 as j varies in 
the range < j < 0.996. The top panels refer to the SN formalism, the bottom panels to the Teukolsky formalism. Linestyles 
are the same as in Fig. |3l 



eigenvalue Aim in ()4.11|) using Leaver's method fn'l. The cigenfunction X,^°^(rf ) can be obtained combining Eq. (IBSp 
and (jC3p . The values of i3„ and ^out can be read off Tables IIIII and IIVI respectively. For increased accuracy we 
evaluate the 7/m's numerically from their definition (|3.1ip . integrating the (complex) angular QNM eigenfunctions 
(see Appendix IB ip . but the numerical results differ only marginally from the simple analytic approximation (|3.12p . 

For fixed {I, m) and n ~ the black hole's response has a relative maximum when the initial data are located very 
close to the horizon. As increases the wave amplitude decreases monotonically, reaching a minimum at ~ 0.86 
(the exact location is only weakly sensitive to the particular mode we consider). This radial location can be seen as 
corresponding to a "node" of the QNM eigenfunction. We also show wave amplitudes for initial data located relatively 
far from the horizon, at = 2. These amplitudes are already quite close to the amplitudes for initial data "falling 
from infinity" . 

The association of nodes with the eigenfunctions is more problematic and ambiguous for larger damping. For 
example, if we consider the fundamental mode of a Schwarzschild black hole (for which Re(a;„) and Im(ti;„) are 
comparable) the wave amplitude turns out to be only weakly sensitive to the location of initial data. The same 
happens when we consider counterrotating modes of (slowly or rapidly rotating) Kerr black holes. 

A natural next step is to study non-localized initial data. We consider static initial data given by a gaussian 
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FIG. 8: Modulus of the SN QNEFs for the fundamental mode and first overtone with I = m = 2. Solid (black) lines are for 
gravitational perturbations, dashed (red) lines for scalar perturbations. 




FIG. 9: Time-domain single-mode waveforms for scalar perturbations of a Kerr black hole with j — 0.98. Each panel shows 
results for diff'erent values of (l,m,n). Dotted (red) lines correspond to initial data located extremely close to the horizon, at 
= 0.6. Solid (black) lines are obtained for initial data localized relatively far away from the black hole, at = 2. Dashed 
(green) lines show the response for initial data localized at the point of minimal excitation, that we determined to be: ~ 0.86 
(Z = m = 2, n = 0); ~ 0.87 (Z = m = 2, n = 1). 



wavepacket 

X(")(to,r) =fce-^('-*-'-^)', (4.14) 

where k and 6 are arbitrary constants. For Schwarzschild black holes, the scattering of gaussian initial data was first 
studied numerically in the time domain by Vishveshwara [lOj. An interpretation of the results in terms of QNM 
expansions for initial data in the far zone was later provided by Andersson [17]. Following [17.] , for initial data in the 
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far zone (large rf ) we get 



-Re 



Re 



n=0 * 



■iu„(t-tBtapt) 



If the gaussian is centered very close to the black hole (rf ^ 0) we have 

i{u!n — 2mfl)kB. 



-Re 



s /TT 

5^" 



"21 

— ^—iuJn{t-tsta.Tt) 



(4.15) 



(4.16) 



Once again, frame dragging effects show up only for gaussians located near the horizon. The absolute value of each 
term in the sum is proportional to exp[— Re(w)/46], with w = cj,^ or w = {un — mfi)^ depending on the location of 
the initial data. Minimizing with respect to h we see that the maximal excitation corresponds to 6 = Re(w)/2: we 
found a similar result in the vibrating string example [recall the discussion below Eq. (|2.22p ]. 

In the limit of narrow gaussians we expect to recover the results for localized initial data. Indeed, if we use the 
(5-function representation setting k = {2y/¥e)^^ and h = l/(4e): 



<5(r» 



1 _(ii^ii)f 



(4.17) 



by taking the limit e ^ we recover (|4.12l) and (|4.13p . 



D. Effective amplitude and initial data 



The present calculation of the QNEFs is only a first step towards the determination of the QNM content of a 
waveform. The degree to which a given QNM is excited depends crucially, through the integral appearing in (13.91) . on 
the initial data /(w,r). Can we identify some physical quantity which is only weakly dependent on the initial data, 
and therefore useful to discuss the detectability of QNMs? 

Andersson and Glampedakis proposed to use a QNM "effective amplitude" defined as follows [H, Con- 
sider a typical ringdown waveform ^. Matched filtering increases the gravitational wave amplitude by VlV, 
N (X Re(ijj„)/Im(ijj„) being the muiiber of cycles, yielding an effective amplitude ^E'cff ~ VN'^. For a single mode 
^' oc C„, so that '^cS ^ VnCu- Now, if we happen to know that the convolution of the initial data with the 
homogeneous solution is reasonably independent of Un, we could write 

*off - VM^n)/lm{uJn)Bn , (4.18) 

which would be the sought criterion for detectability of a given mode. In [l^, [l^] the real part of (|4.18p was referred 
to as the "effective gravitational wave amplitude" , and used to discuss QNM detectability. 

The derivation of Eq. (|4.18p relies upon the assumption that the proportionality factor between the C„'s and the 
BnS is weakly dependent on aj„, for only in this case is the previous expression valid. Eqs. (I4.12p and (|4.13p show 
that such an assumption is reasonable for localized initial data. Unfortunately this won't be true in more general 
situations. A simple counterexample is provided by gaussian initial data. According to (|4.18p QNM detectability 
should be weakly dependent on cj„ (at least for low- lying modes), but we know for a fact that the correct answer 
for gaussian initial data, Eq. (|4.15p . has an exponential dependence on tj„. Any notion of detectability based on the 
naive argument leading to (j4.18p is therefore dangerous, and a case-by-case analysis of the initial data is mandatory. 
A study of the black hole's response to physically reasonable initial data will be the subject of future work. 



V. CONCLUSIONS AND OUTLOOK 



Motivated by the prospect to test the Kerr nature of astrophysical black holes by gravitational wave observations, 
in this paper we compute QNEFs for general-spin perturbations of Kerr black holes. For corotating modes with I = m 
we find that QNEFs tend to zero in the extremal limit, and that the overtone contribution is more significant when 
the black hole is fast rotating. This result is confirmed by numerical time-evolutions of scalar perturbations of Kerr 
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black holes [i^l ■ In Appendix [X] we also present the first analytical calculation of the large-n QNEFs for static black 
holes, including the Schwarzschild and Reissner-Nordstrom metrics. 

QNEFs are "universal" properties of the Kerr metric: they do not depend on the physical nature of the perturbation. 
Suppose that a distorted Kerr black hole is formed as a result of a binary merger. The detailed physics of the merger, 
depending on parameters such as the masses, spin magnitude and inclination of the binary members, will affect the 
relative overtone excitation through the initial data function /(cj , r) appearing in the definition ()3.9p of the excitation 
coefficients. Therefore the problem of determining the relative overtone excitation in ringdown waveforms reduces to 
the determination of the function I(uj , r) from numerical simulations. 

Of course the present formalims relies on the possibility to use linear perturbation theory to extract ringdown 
waveforms from numerical simulations. In the nonlinear regime the total mass and angular momentum of the dynam- 
ical spacetime could deviate significantly from their values for the final Kerr black hole, introducing systematic (and 
perhaps time-dependent) redshifts in the QNM spectrum [4l'|. An invariant criterion to monitor the applicability 
of linear perturbation theory is to use the Petrov classification, looking at the relative deviation of the background 
spacetime from Type D. The Petrov type of the background can be determined, for example, using a wave-extraction 
formalism based on the quasi-Kinnersley tetrad and looking at deviations of the scalar invariant S from the value 
(5 = 1) that it would have for a Kerr black hole \A2^ . 

A significant technical challenge in the calculation of the excitation coefhcients C„ is the renormalization of the 
divergent integral p.9p for generic initial data. We will address this problem in future publications. As a tcstbcd 
we plan to study simple models amenable to perturbation theory (eg. the old problem of particles plunging into 
Kerr black holes). Then we wish to explore the nonlinear regime in relatively symmetric situations, such as head-on 
collisions or nonlinear simulations of single, distorted black holes. Finally we will consider the harder, astrophysically 
realistic problem of determining the relative overtone excitation from nonlinear simulations of merging black holes, 
using Newman-Penrose based wave extraction of the ringdown waveform from full numerical relativity simulations. 

An accurate determination of the relative QNM excitation in realistic binary mergers would have interesting as- 
trophysical applications. For example it could be used to assess the feasibility of tests of the no-hair theorem with 
Earth-based and space-based gravitational wave interferometers, and to improve present estimates of the "ringdown 
braking" phase observed in numerical and analytical calculations of gravitational wave recoil [H, . 
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APPENDIX A: ASYMPTOTIC QUASINORMAL EXCITATION FACTORS OF STATIC BLACK HOLES 

One of the fundamental prerequisites in order to have a well-defined QNM expansion is that the QNEFs -B„ 
"converge" in some sense. So far this problem has not received much attention. A numerical study of the scalar 
QNEFs for large n and non-rotating black holes can be found in p^ . In this Appendix we study analytically the 
large-n behavior of the QNEFs for Schwarzschild black holes and other static geometries. We explicitly compute the 
dependence of the QNEFs on n (including subleading effects) in three important cases, including the Schwarzschild 
and Reissner-Nordstrom spacetimes and the widely used Poschl- Teller approximation of the Regge- Wheeler potential. 
We also show that for a general class of static spacetimes the leading behavior of the QNEFs for large n is ^ 1/n, and 
we conjecture that the same asymptotic behavior should apply to Kerr black holes. It should be easy to generalize 
our technique to compute subleading effects. Finally, as a simple counterexample to the generality of this behavior, 
we show that the QNEFs for a potential barrier scale as for large n. 

1. Schwarzschild black holes 

Our general method relies on the monodromy argument by Motl and Neitzke p5| . devised to compute highly 
damped (i.e. large-n) QNM frequencies (see [S^ and for reviews). Using this method Neitzke computed the 
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reflection and transmission coefficients R and T for large n [47|. From his calculation we can simply read off the 
in-going and out-going wave amplitudes for the Regge- Wheeler equation: 

_ 1 _ e^"-H + 2cos(^g) 

^ T ~ e/3c^ _ 1 ' ^^^> 

where in our units (G = c = 2M = 1) /3 = 47r is the inverse of the black hole's Hawking temperature. This result 
is valid for large n and for any angular quantum number I. QNM frequencies, being poles of the reflection and 
transmission coefficients, satisfy 

e'^"" + l + 2cos(7rs) =0. (A3) 

The calculation of the QNEFs Bn for large n is now straightforward. Denoting the large-n limit of the QNEFs for 
spin-s perturbations by Bn'' we find 

(.) ^ / ^out \ ^ »cos(7rs/2) 

" - {2ujd^AiJ ^^^^ piOr^eP-- ■ ^ ' 

This leads to the result anticipated in the main text, Eq. (|3.17p : 

b(°) = = L ^ (n ^ oo) . (A5) 

3[±log3- (2n-hl)7ri] ' ^ ' ^ ' 

This is our main result for the Schwarzschild QNEFs in the large damping limit. The analytic prediction agrees very 
well with Figure 2 in T^, where the QNEFs are plotted for the first ~ 200 scalar modes. The calculation of Ref. [T8| 
shows that for large n QNEFs are independent of I and drop linearly with n, in agreement with our prediction. Notice 
also that for electromagnetic perturbations the asymptotic QNEFs Bn would be zero: higher-order corrections 
must be taken into account to obtain a nontrivial result. 



2. Reissner-Nordstrom black holes 

The monodromy technique is easily adapted to other black hole backgrounds [i^ , allowing the computation of the 
asymptotic QNEFs for a wide class of geometries. Take for instance the Reissner-Nordstrom metric, describing a 
black hole with mass M and charge Q. We follow Motl and Neitzke [i^ and fix units such that M = (fc -|- l)/2 and 
= k. With this choice the outer horizon is located at r = 1, and we get [13] 

A - ^- + 2 + 3e-^-" 

^ T~ eP'^ - 1 ' ^ ^ 

Aout = ^^±*^/3^t__, (A7) 

for any massless integer-spin field (at odds with the Schwarzschild case, where this symmetry is broken for electro- 
magnetic perturbations). Here /3 = 47r/(l — fc) and j3i = —k^fi denote the Hawking temperatures of the outer and 
inner horizons, respectively. For large n the QNM frequencies a;„ are solutions of 

e'''^ + 2 + ie-'^"^ = . (A8) 

An elementary calculation similar to the Schwarzschild case yields 

Blf> = 3 3 , (n-^oo). (A9) 

In (jA7[) and (jA9p the plus sign refers to electromagnetic-gravitational perturbations, and the minus sign refers to 
scalar perturbations. The fact that in the large-n limit the predictions for charged black holes do not reduce to the 
Schwarzschild predictions as Q — *■ is known, and agrees with numerical results [s^. In general, the transcendental 
equation (|A8|) must be solved numerically to obtain the QNM frequencies. As in the Schwarzschild case, one can 
verify that the QNEFs B„ ~ 1/n. 
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3. The Poschl- Teller potential 



The Poschl- Teller potential has been used in many studies of black hole perturbation theory as a useful approxi- 
mation of the Regge- Wheeler potential. The idea is to replace the "true" potential barrier by an analytic expression 



of the form V 
closed form |3t 



V 



PT 



Vq cosh a{x — Xq), which has the advantage that many calculations can be carried out in 



|49| . From the analytic expression of the reflection and transmission coefficients given in [36[ we get 



^out 



r(l + iuja) r(— iw/a) 

r(i + /3)r(-/?) ' 

r(l -hiw/a)r(iw/a) 
r(l + /? iu;/a)r(-/3 + icj/a) 



(AlO) 
(All) 



where /? = — 1/2 + ^^1/4 — Vq/o^. As usual, the QNM frequencies are poles of Ai^. They can be obtained imposing 
(say) r(— /3 + icj/a)^^ — (or equivalently — /3 + iuj/a = —n). Since i3„ ^ (^i^^in)"'^, the only nonzero contribu- 
tion to Bn must come from the term [r(l -I- iuo/a) r(ia;/Q;)/r(l + (3 + iu/a)] 9i^(l/r(— /3 -I- iiu/a)). Introducing the 
polygamma function V'Oi such that r'(a;) = T(x)iPq{x), the derivative we need is 



d 

dco 



it/jo{—n) (— l)"in! 
Q;r(— n) a 



and the final result for the asymptotic QNEFs of the Poschl- Teller potential reads 

m(-l)"+i r{n- f3)T{l + 2f3-n) 



B. 



PT 



2L^„r(i + /3)r(-/3) 



nW{(3 



(A12) 



(A13) 



Once again, asymptotically B 



PT 



1/n. 



4. Generic static black hole spacetimes and other possible extensions 



An investigation of generic spacetimes with the present technique should be carried out case by case [48i |. However, 
if we are only interested in the leading-order behavior for l arg e n we can content ourselves with a simple Born 
approximation to the scattering amplitude (see for instance [50|) to deduce that cu ~ —in, so that generically the 
QNEFs B„ - 1/n for large n. 

The generalization of the present results to Kerr black holes is more difficult. Previous numerical calculations of 
the asymptotic QNM frequencies [H, show that the imaginary part is still proportional to n, so we expect that 
(to leading order) the Kerr QNEFs should still scale as B^°" ^ 1/n. 



5. A potential barrier 

As a last example we consider the QNEFs for a potential barrier of height Vq and width x, i.e. the wave equation 
F(a;)]* = , with potential 

V{x) = I ^f" < n ^ - (A14) 

^ ' 10 for X < and x > x . ^ ' 



Defining k — y^ui'^ — Vq, the solution in each of the three different regions is 

* = 6-^'^"= , X <0, (A15) 

* = Ae*'==" -I- Be"''"'' , 0<x<x, (A16) 
^ ^ (jgii^x ^ jjg-ii^x ^ 0<a;<x. (A17) 

The constants A, B, C, D are obtained by imposing continuity of the field and its derivative at a; = and x — x, so 
that A + B — 1 and —oj = kjA— B)/ (A + B). QNM frequencies are such that D{u!) = 0. This problem was considered 



by Chandrasekhar and Detweiler [5l|, who showed that for large mode number n the QNM frequencies behave as 



(7i7r — iuji), with tUi a solution to the transcendental equation ~ —2uJie^^/^. We have solved the differential 
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equation numerically, and our results agree very well with this prediction. For large overtone n we find the analytical 
result 

C ~ A.e-"('^+'=) (1 - e^"'^) , (A18) 

We also get 

At the QNM frequencies VJj^e^"'^ - 16w^ = 0, and we get the following QNEFs 

^barrier ^ _C Vo .^21) 

where we used the fact that in this regime Im(a;) — oo and e~"('^+'^) ^ 0. The result (|A21[) is in good agreement 
with numerical results and shows that the 1/n behavior, although very general, is not universal, for in this case the 
QNEFs scale as 



APPENDIX B: DETAILS ON THE CALCULATION OF THE QUASINORMAL EXCITATION FACTORS 

This Appendix provides details on the calculation of QNEFs for general-spin perturbations of Kerr black holes. 
Most of the material can be found in the original papers by Leaver [ll|, [H, (see also [1, HI, HJ), but we find it 
convenient to summarize here the equations which are necessary to implement the computational procedure. 



1. Expansion of the angular wavefunction 

If we express the angular wavefunction Sim as 

oo 

Sim{u) = e-^^il + (1 - u)^ + ' (Bl) 

n=0 

the coefficients must satisfy the following recurrence relation 

a'oal + pyo = 0, (B2) 

a'a'+i+/?X+7'4-i = 0, (B3) 

where 

ai = -2(n + l)(n + 2^1 + 1), (B4) 
13^ = n{n -l) + 2n{ki + fcz + 1 - 2aLo) - 

[2auj{2ki + s + 1) - (fci + k2)iki + fca + 1)] - a^ui^ - s{s + 1) - A,™ , (B5) 

7^ = 2aw(n + fci + fca + s) . (B6) 

Regularity of the solution at the boundaries implies that the sequence of expansion coefficients must be minimal, and 
that the separation constant Aim must be a root of the continued fraction 



_ .e «57? a!72 «27l (^^. 



or any of its inversions. 
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2. Jaffe expansion of the radial wavefunction 

As discussed below Eq. (|4.9[) . to compute the Kerr QNEFs we need an accurate representation of the solutions of 
the radial equation. A convenient series solution close to the horizon can be found by methods due to Jaffe (see 

OG y K n 

n=0 ^ 

where the notation is the same as in Sec. IIV A1 The coefficients aj^ are normalized so that — 1, consistently with 
Eq. (|4.7ap . They can be obtained from the recurrence relation 

« + f3o< = 0- (B9) 
a>;+i+/3X+7;a;,i = 0, (BIO) 



where 



and 



a; = + {co + l)n + co, (Bll) 
/3; = -2n2 + (ci+2)n + C3, (B12) 
Yn - n2 + (c2-3)n + C4-C2 + 2, (B13) 



Co — 1 — s — iuj (^— — arnj , (B14) 

ci = -4: + 2iuj{2 + b) + j(^^-am^ , (B15) 

C2 — s + 3 — 3iuj — — (^— — arri^ , (B16) 

C3 = uj'^{i + 2b- a'^) - 2amuj - s - 1 + {2 + b)iuj ~ Ai„, + ^^^^^{'^ - am) , (B17) 

C4 = 5 + 1 - 2cj^ - (2s + 3)zcj ~ ^ ~ Q'^) ■ (B18) 



QNM frequencies satisfy a continued-fraction relation analogous to (jB7p . with the superscript replaced by a super- 
script r. Given the QNM frequencies and the corresponding angular eigenvalues Aim it is a simple matter to compute 
i?r+ for any finite, not very large value of r. Unfortunately, the convergence of the Jaffe expansion gets worse for 
large values of r. This is precisely the region where the wavefunction can be expressed as a sum of in- and out-going 
components, a necessary procedure to evaluate Af^ and A^^^. and their derivatives. For this reason we must resort to 
the Coulomb wavefunction representation first introduced by Leaver [T^ ITsj , and summarized below. 

3. Coulomb wavefunction expansion of the radial wavefunction 

Let us first introduce a new wavefunction h, related to the Teukolsky wavefunction R by 

R={r- r_)-i+*^+ (r - r+y-"'^ h{z) . (B19) 
where z = aj(r — r_). The wavefunction h satisfies the generalized spheroidal wave equation 

z{z - ujxo) [d^ h+{l- 2r]/z)h] + CiLud^^h + (C2 -I- Cslo/z) h^Q, (B20) 

with 

xq — h , 1] ~ —is — to , (B21) 

Ci = b{l - s - 2ic7+) , (B22) 

C2 = -Aim - s{s + 1) + ibsLU + 2u? - c?u? + bup- , (B23) 

C3 = 6(-l + s + iw)(l + iw-2i(T+). (B24) 
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We next expand h as 

h{z)= aLUL+uiv,z), (B25) 

L — — C30 

where Ur, 4^,^(11, z) is any combination of the special functions F]^^^(r] , z) and GL-i-^{ri , z), known as Coulomb wave- 
functions (see (l^ . [5^ for more details on these functions). We are interested in two particular combinations of 
Coulomb wavefunctions: 

U+^^{r^,z) = GL+.{Tl,z) + iFL+,{ri,z), (B26) 
^L+.(^>^) = GL+.{r],z)~iFL+,{ri,z). (B27) 
As z ^ 00 these functions have a simple asymptotic behavior j52l | : 



(B28) 



with 



Other properties of these functions are given in Appendix IB 41 The coefficients in (|B25p are determined by 

aiOi+i + /^lCl + 7LaL-i = , (B30) 

where 

[L + 1 + v)(2L + 2v + i) 

/3l = {L + ,y){L + iy+l) + C2+ i^L + i,){L + v+\) + ^)(^ + ^ + " ^1 ^ ^s] , (B32) 



7. = - ^Lf^^tl lt + -)(^ + - - 1)^0 - Cr(L + . - 1) - C3] • (B33) 



It is convenient to define the quantities [52 



From these definitions it follows that 



Rl^^, Ll^^. (B34) 



, Ll = --—^. . (B35) 

For a minimal solution of the recurrence relation the (as yet undetermined) parameter v must be a root of [T^ 

a-i7o a-27-i "-37-2 , ao7i "172 "273 /-Doc^ 
'^" = ^ ^ R ■■■ + ^5 — ^5 — 1^ — ■■■■ ^^^^> 

P-l- P-2- P-3- Pi- P2- P3- 

The parameter ly can also be determined as a solution of the equation jssj 

RlLl-i = 1. (B37) 

However, as suggested in (53 |. the most convenient numerical choice is to solve for 

pL+aLRL+i+lLLL-i = 0. (B38) 

We are interested only in the solutions that map to the correct asymptotic value v = I as lo (corrections of order 
can be found in [ssj). Roots z/ that are integer multiples of 1/2 are usually spurious and must be discarded 1^. 
For slowly damped modes, physically meaningful roots are most easily obtained setting L = OorL = — lin (|B38[) . 
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We verified that different numerical procedures yield excellent agreement on the resulting values of v. For reference, 
we list some of these values along with the corresponding QNM frequencies (for a = 0) in Table rvTIl 

The procedure to compute the functions i?oo+ and i?oo_ is now straightforward, at least in principle. Compute the 
parameter v solving (|B36p . (jB37p or (|B38p . Then compute the coefficients for positive and negative L from (jB34p . 
setting for instance oq = 1 (the final results arc of course unaffected by this or other normalization choices). Finally, 
the general solution follows from (|B19p and (|B25p . Imposing the normalization condition ()4.8bp we get 

i?oo±=^ ^ ^ E aLUt^,, (B39) 

L — — oo 

where from (jB28p the normalization constants read 

oo 

K±— ^ Q^g±[-«';iog(2i^)-«(i'+'^)^/2+i<Ti] ^ (B40) 

L — — oo 

The nontrivial step in the calculation of (|B39p is the evaluation of the Coulomb wavefunctions Uj^^^ for complex 
arguments. Computing these functions is particularly tricky for large values of \L\, which must be included for (jB39P 
to converge with acceptable accuracy. Typically a good precision (of the order of five or more significant digits) is 
achieved summing terms up to \L\ ~ 15. 

We verified our results in a number of ways. We first computed the functions U^^^{rj, z) using their relation with the 
confluent hypergeometric function U, Eq. (125) in jl2i] . Then we checked their values using a similar representation in 
terms of the confluent hypergeometric function M, that can be derived from the integral representation (111) of p^ . 
The reason for switching to U and M is that confluent hypergeometric functions (unlike Coulomb wavefunctions) are 
implemented in the present version of Mathematica. Our Mathematica calculations agree very well with a Fortran 
subroutine to compute Coulomb wavefunctions for complex arguments [H^, that was eventually used to obtain the 
results in this paper. Further checks on our numerical results are discussed in the following section. 



4. Numerical checks of the solutions: some useful identities 



Here we describe a number of consistency checks we performed on our numerical solutions. We performed the 
calculations using both Mathematica and a Fortran code. Our final results on the QNEFs typically agree to one part 
in 10^ or better. 

The convergence of the series expansion (|B8p was checked via a direct high-precision integration of the Teukolsky 
equation. When evaluated in the matching region (typically we use r 5 — 7) the two methods agree to one part in 
10^ or better. 

Besides checking the validity of the asymptotic expansion (|B28I) , the following equalities proved useful to determine 
the accuracy of the functions Ui^'^'^- 



d'U^^^ , / 277 {L + u){L + ly + 1) 



= (B41) 



^(i + VY + T]^ 

(L + .) v/(lT7TTFT^^±+„+, = (2(i + .) + 1) + iL + .){L + .+ l) ^ 
{L + u + iy{L + vY+i^^Ut^^^^ (B43) 

(i + ^ = v/(rWT^^f+._, - { + ^] Ut^^ (B45) 

(B46) 



Numerical calculations of U^^'^ usually fail for large \L\. This is true both when we use their representation in terms 
of confluent hypergeometrics and when we use the Fortran routines of Ref . [55| . In the sum (jB39p we included only 
those values of L for which the identity (|B42p is satisfied (in modulus) to better than one part in 10"*. 
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An important test on the hnear combinations Roo± is based on the calculation of their Wronskian. One can prove 
easily that the Wronskian W between any two solutions of the Teukolsky equation satisfies 



W{r) = iy(ro)exp 



{x - r+){x - r^) 



(B47) 



where vq is any fixed point. Considering the two solutions i?oo+ and -Roo_ we can evaluate the constant W{ro) by 
computing the Wronskian at infinity from their asymptotic behavior (j4.8bp . with the result 

W{r„) = 2tuj [(ro - r+)(ro - r_)]-^'+'^ . (B48) 

Our wavefunctions satisfy this relation at different selected values of rg. A final test consists in evaluating A^ut by a 
direct integration of the Teukolsky equation and comparing the results with the method described in the main text: 
the agreement is at the level of one percent or better. 

APPENDIX C: TRANSFORMATION BETWEEN THE TEUKOLSKY AND SASAKI-NAKAMURA 

WAVEFUNCTIONS 

In the main text we computed the QNEFs using the Teukolsky formalism. Quite often it is com put ationally 
convenient to use the equivalent formalism developed for s = —2 by Sasaki and Nakamura p^ . |23|. |24'| and its 
generalization for other spins [25]. In general, switching between the two formalisms is not trivial. Fortunately, 
deriving the asymptotic behavior of the respective wavefunctions and the relation between the QNEFs is quite simple, 
as we show below. 

We follow Sasaki and Tagoshi [23| and denote the SN function by X'-'^K At infinity, the asymptotic behavior of 
X'-^^ is the same for all s: 

lim X(^) ~ Aine-^'^-* + Aoute'"''* , (CI) 

r — >oo 

The normalization at the horizon is fixed by our choice for the Teukolsky function, Eqs. (|4.7aP and (|B8p . We will 
define the QNEFs in the SN formalism in the usual way: 



2uj V cLj 



(C2) 



Below we establish the relation between B^^^ and the QNEFs b!^'' in the Teukolsky formalism, as computed in the 
main text. 

1. Scalar perturbations 

In the scalar case, the generalized SN function is related to the Teukolsky radial function R by [2^ 

X(") = (r2+a2)i/2i?, (C3) 



and satisfies 



where 



^X(")-Z^(")X(0)=0, (C4) 
dri 



^(0) ^ _ K^-iA,n ~ 2arn. + a^.^)A ^ ^ ^ (C5) 

(r^ + a'^Y ar^ 



with K ^ {r^ + a?)uj ~ am, and G — rA{r'^ + a?) ^. 
From Eq. (jC3p we have immediately 



- . (C6) 
The normalization of the Teukolsky function, Eq. (j4.7ap . implies that at the horizon 

lim X^"^ ^ {rl+a'^)^/'^ir+ -r-)-^+"^+"'+e"^'''+{r -r+)-"'+ . (C7) 

r— s-r+ 
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2. Electromagnetic perturbations 

For the electromagnetic case, the generahzed SN transformation is 

= ^J^^ [a + Pdr R^-'^) , (C8) 

where we can choose [1^ 

From the asymptotic behavior of the Teukolsky function 

Hm - Alr-^e-''^'-' (1 + ^n/r) + tre*'"''* (1 + M2/r) , (CIO) 

r — ^oo 

with ^1 ~ [—2uj + i{2 — Aim — a^w^)] / {2uj) and /12 = [« (^/jn + a-^w^)] / {2uj) we get 

^in = ~T7~^in; ^out = ^ 7; 7 2"^"^°"' ' (CH) 

2^a; Zamuj — Aim — a'^Lo 



and finally 



4-1) ^ ^B(-i) . (C12) 



Similarly, from the asymptotic behavior of the Teukolsky function at the horizon (|4.7ap we find 

lim X(-i' = (l-2i(T+)rI^^^(r+-r_)*"+*'^+e*""^+(r-r+)-""+. (C13) 

r — >r+ 

3. Gravitational perturbations 

For the gravitational case the transformation between the SN and Teukolsky wavefunctions is more complex, but 
here we only need the asymptotic behavior (for more details we refer to the original papers (2^ . [23j and to the review 
in d^l)- One can show, eg. from Eqs. (65)-(69) in that asymptotically the Teukolsky and SN amplitudes satisfy 

" "4^^'" ' " " A(A + 2) - 6zw - 12aw(ac^ ~ m) ' 

so that 

4"'^ - T7^ ^ 1 -.B^-'^ • (C15) 

A(A + 2) — QiLO — 12auj[auj — m) 

with A = Aim + {a.'^y' — 2amuj. Expressing Eqs. (65) and (68) in ^24.] in Leaver's units and combining them with 
Eq. (|4.7ap . the normalization at the horizon is: 

lim ^ d{r+ - ,._)-i+i'^+»'^+e^"'^+(r - r , (C16) 

r — 

where d = [(8 — 12ia; — Auj'^)r\ + {12iam — 8 + Samuj + 6icj)r+ — Aa^ni? — 6iam + 2] . 
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TABLE II: Angular separation constants Ai„ 
parameter and for some values of m. 



of the first two overtones with s = and Z = 2, for several values of rotation 



0, n = 



J 



-1 



: 



0.00 6.0000+O.OOOOOi 6.0000+O.OOOOOi 6.0000+O.OOOOOi 6.0000+O.OOOOOi 6.0000+O.OOOOOi 

0.20 5.9989+0.00050i 5.9964+0.00156i 5.9953+0.00196i 5.9959+0.00166i 5.9985+0.00057i 

0.40 5.9959+0.001901 5.9862+0.00604i 5.9807+0.00781i 5.9819+0.00682i 5.9930+0.00243i 

0.50 5.9939+0.00289i 5.9788+0.00929i 5.9695+0.01217i 5.9702+0.01080i 5.9880+0.00392i 

0.60 5.9917+0.00405/ 5.9700+0.01317i 5.9553+0.01743i 5.9544+0.01570'j 5.9808+0.00582't 

0.80 5.9865+0.00684i 5.9482+0.02266i 5.9167+0.03030i 5.9056+0.02776't 5.9547+0.01063i 

0.90 5.9836+0.00845i 5.9352+0.02817i 5.8911+0.03744i 5.8674+0.03351i 5.9294+0.01271j 

0.98 5.9812+0.00982i 5.9237+0.03289i 5.8668+0.04293i 5.8245+0.03416i 5.8884+0.01030i 



0, n : 



m = -2 



m = —1 



m = 



m = 1 



m = 2 



0.00 6.0000+0. 

0.20 5.9994+0. 
0.40 5.9982+0. 
0.50 5.9974+0. 
0.60 5.9967+0. 
0.80 5.9954+0. 
0.90 5.9948+0. 
0.98 5.9944+0. 



OOOOOi 

00147i 
005521 
00836i 
01169i 
01961i 
02414i 
02799i 



6.0000+0. 

5.9980+0. 
5.9927+0. 
5.9890+0. 
5.9845+0. 
5.9732+0. 
5.9662+0. 
5.9599+0. 



OOOOOi 

00455i 
017621 
02709i 
03838i 
06591i 
08181i 
09535i 



6.0000+0. 

5.9973+0. 
5.9886+0. 
5.9815+0. 
5.9720+0. 
5.9439+0. 
5.9237+0. 
5.9043+0. 



OOOOOz 

00575?: 
02294?: 
03577i 
051294 
08920i 
10993i 
12559i 



6.0000+0. 

5.9975+0, 
5.9882+0, 
5.9797+0, 
5.9676+0, 
5.9256+0, 
5.8894+0, 
5.8509+0, 



OOOOOi 

00486?: 
02010? 
03186't 
04641? 
08210? 
09850i 
09778i 



6.0000+0, 

5.9991+0, 
5.9951+0, 
5.9911+0, 
5.9850+0, 
5.9606+0, 
5.9347+0, 
5.8908+0, 



OOOOOi 

00168?: 
00719? 
01163't 
01727? 
03168? 
03791i 
03067i 



TABLE III: QNEFs B^°^ of the first two overtones with s = and I = 2, for several values of rotation parameter and for some 
values of m. 

s = 0, n = 



J 



m — —1 



m — 



m — 1 



m = 2 



0.00 0.11936+0.01343i 0.11936+0.01343i 0.11936+0.01343i 

0.20 0.12028-0.00281i 0.12013+0.00695z 0.11907+0.01746i 

0.40 0.12079-0.01047? 0.12090+0.00828? 0.11714+0.03012? 

0.50 0.12156-0.01146't 0. 12130+0. 01197i 0.11447+0.04016i 

0.60 0.12283-0.01049't 0. 12149+0. 01794't 0.10945+0.05296i 

0.80 0.12720-0.00092? 0.11939+0.03899? 0.08459+0.08650i 

0.90 0.13027+0.01075i 0.11456+0.05709i 0.05691+0.10467i 

0.98 0.13253+0.03151i 0.10313+0.08251i 0.01644+0.11542i 



0. 11936+0.01343i 0. 11936+0.01343i 

0. 11687+0.02868i 0. 11332+0.04050i 

0.10769+0.05442? 0.09045+0.07951? 

0.09715+0.07141't 0.06518+0.10076? 

0.07919+0.09011? 0.02441+0.11621? 

0.00372+0.11379i -0.09556+0.05169i 

-0.05609+0.08585i -0.05739-0.07009i 

-0.07590+0.01089i 0.03015+0.03247i 



s ■■ 



0, n = 1 



m = -2 



m = —1 



m = 



m = 1 



m = 2 



0.00 0.03552-0.26427i 0.03552-0.26427i 0.03552-0.26427i 0.03552-0.26427i 0.03552-0.26427i 

0.20 0.00031-0.23879? 0.02062-0.25361? 0.04536-0.26760? 0.07508-0.27998? 0.11026-0.28969? 

0.40 -0.01305-0.22194i 0.02373-0. 25063i 0.07830-0.27535i 0.15482-0.28823? 0.25502-0.27656? 

0.50 -0.01418-0.21735? 0.03186-0. 25195i 0.10651-0.27817? 0.21694-0.27760? 0.36021-0. 21856i 

0.60 -0.01204-0.21552? 0.04520-0.25482? 0.14519-0.27668? 0.29610-0.24132? 0.46944-0.08037? 

0.80 0.00433-0.22175Z 0.09487-0.26260i 0.26097-0.23285i 0.44182+0.01597i 0.24418+0.54878i 

0.90 0.02430-0.23258Z 0.14156-0.26409i 0.33334-0.16103i 0.33715+0.26654i -0.50193+0.31383z 

0.98 0.06291-0.25169i 0.21530-0.25651i 0.39147-0.04418i 0.06435+0.29230i 0.23752-0. 12729i 
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TABLE IV: Amplitudes ^out of the first two overtones with s = and I = 2, for several values of rotation parameter and for 
some values of {I ,m). 

s = 0, n = 



J 



-1 



m ■ 



: 



m • 



1 



m • 



0.00 
0.20 
0.40 
0.50 
0.60 
0.80 
0.90 
0.98 



1.1472-1.1426i 
1. 1068-1. 2514i 
1.0516-1.4122i 
1.0057-1.5269i 
0.9311-1.68021 
0.5060-2.2157i 
-0.3226-2.7081i 
-3.8026-2.0011Z 



1. 1472-1. 1426i 
1.1252-1.2067i 
1.0888-1. 3208i 
1.0547-1.4093't 
0.9955-1.5331i 
0.6304-1.9952i 
-0.1229-2.4419i 
-3.3378-1.6360i 



1. 1472-1. 1426i 
1. 1455-1. 160H 
1.1353-1.22111 
1. 1203-1. 2776i 
1.0886-1.3636i 
0.8437-1.7289i 
0.2580-2.1461i 
-2.6005-1.6232i 



1.1472 
1.1677- 

1.1904- 
1.2013- 
1.2085- 
1.1526- 
0.8860- 
-1.1381 



1.1426i 
1.1115i 
1.1119i 
1.1287i 
1.1641j 
1.3741i 
1.7127i 
■2.0428i 



1.1472- 
1.1915- 

1.2522- 
1.2930- 
1.3444- 
1.5014- 
1.6345- 
1.7179- 



1.1426i 
1.0608Z 
0.9921j 
0.9595j 
0.9272i 
0.8727i 
0.9057i 
1.479H 



j m = —2 
0.00 -2.9403-|-3.1839i 
0.20 -2.3498+2.9311i 
0.40 -1.8312+2.6818i 
0.50 -1.5742+2.5630i 
0.60 -1.3002-|-2.4481i 
0.80 -0.5767-|-2.2022i 
0.90 0.0721-M.9681i 
0.98 1.2103-F0.9694i 



m = —1 
-2.9403-h3.1839i 
-2.6009-F3.0411i 
-2.2132-h2.8695i 
-1.9831+2.77731 
-1.7093-F2.6818i 
-0.8741-F2.4579i 
-0.0536-^2.2109* 
1.3828-|-0.9603i 



I = 0, n = 1 
m = 
-2.9403-h3.1839i 
-2.8990+3. 1573i 
-2.7542+3.0781i 
-2.6206+3.01951 
-2.4198+2.9493i 
-1.5934+2.7666i 
-0.5925+2.5617i 
1.4079+1.1980i 



m = 1 
-2.9403+3.1839i 
-3.2543+3.2790i 
-3.5300+3.2999i 
-3.6326+3.2719i 
-3.6819+3.2125i 
-3.3396+2.99901 
-2.4246+2.8572i 
0.4182+1.9828i 



m = 2 
■2.9403+3.1839i 
■3.6792+3.4045i 
■4.6593+3.5080i 
■5.2745+3.4588i 
■5.9898+3.2826i 
■7.6627+2. 1745i 
■8.3413+0.9505i 
■8.0724+0.3919i 
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TABLE V: QNEFs B^~^^ of the first two overtones with s = —2 and I = 2, for several values of the rotation parameter and all 
values of m. 

Sasaki-Nakamura, s = —2, n = 



m = -2 



m — —1 



rn = 



m = 1 



TO = 2 



0.00 0.12690+0.02032i 0.12690+0.02032i 0.12690+0.02032i 0.12690+0.02032i 

0.20 0.13384+0.00380i 0.13066+0.01398i 0.12644+0.02431i 0.12114+0.03460i 

0.40 0.13951-0.00575?: 0.13409+0.01509?: 0. 12399+0. 03694i 0.10873+0.05810?: 

0.50 0.14255-0.00798? 0.13574+0.01864? 0.12088+0.04702? 0.096848+0.073395? 

0.60 0.14598-0.00838? 0.13716+0.02454i 0.11523+0.05998i 0.077845+0.090076i 

0.80 0.15473-0.00179? 0.13753+0.04581i 0.08743+0.09424? 0.00070+0.10841? 

0.90 0.16055+0.00821? 0.13413+0.06426? 0.05603+0.11183? -0.059511+0.072290? 

0.98 0.16652+0.02731i 0.12487+0.08994? 0.01200+0.11809i -0.056299-0.016701? 



0.12690+0.02032? 
0.11481+0.04465i 

0.089192+0.076810i 
0.066032+0.094125? 

0.03074+0.10817? 
-0.091113+0.061347? 
-0.061104-0.089943? 
-0.042248+0.067331? 



Teukolsky, s = —2, n = 



?n = -2 



m = — 1 



?n = 



m = 1 



?n = 2 



0.00 0.025587-0.016876? 0.025587-0.016876z 0.025587-0.016876? 0.025587-0.016876i 

0.20 0.015858-0.015534i 0.020659-0.016165? 0.026663-0.016183? 0.034089-0.015199? 

0.40 0.010593-0.013107i 0.018158-0.014516i 0.029957-0.013582i 0.047365-0.006739i 

0.50 0.008973-0.011907i 0.017585-0.013461? 0.032435-0.011025i 0.055380+0.003377i 

0.60 0.007817-0.010778? 0.017385-0.012229? 0.035338-0.007035? 0.062135+0.021050? 

0.80 0.006570-0.008740i 0.017957-0.008866i 0.040011+0.008791i 0.034721+0.091829i 

0.90 0.006467-0.007753i 0.018690-0.006243i 0.037266+0.023214i -0.05801 +0.10686i 

0.98 0.006929-0.006758? 0.019607-0.002616? 0.025180+0.038227« -0.10800 -0.041437? 



0.025587-0.016876? 
0. 043137-0. 012654i 
0.070428+0.012901't 
0.081892+0.046677? 

0.07301 +0.10619? 
-0.24081 +0.15010? 

-0.13315 -0.45394i 
-0.60610 +0.26215? 



Sasaki-Nakamura, s = —2, n = 1 



?n = -2 



m = — 1 



?n = 



m = 1 



?n = 2 



0.00 0.04769-0.22379? 0.04769-0.22379? 0.04769-0.22379? 0.04769-0.22379? 

0.20 0. 01275-0. 21222i 0.03286-0.22075't 0.05588-0. 22700i 0.08149-0.23024't 

0.40 -0.00572-0.20420't 0.03229-0.22415't 0.08389-0. 23496i 0.14621-0.22961't 

0.50 -0.01064-0.20281? 0.03752-0.22859? 0.10852-0.23864? 0.19604-0.21680? 

0.60 -0.01302-0.20359i 0.04716-0.23482? 0.14339-0.23886? 0.25927-0.18219? 

0.80 -0.00904-0.21482? 0.08713-0.25154? 0.25286-0.19967? 0.36945+0.05388? 

0.90 0.00093-0.22925? 0.12656-0.26018? 0.32055-0.12532? 0.23328+0.28598? 

0.98 0.02390-0.25495? 0.18872-0.26443? 0.35732-0.00915? -0.08578+0.18065? 



0.04769-0.22379? 

0.10888-0.22995? 

0.20947-0.20386? 

0.27571-0.15875? 

0.34491-0.07025? 

0.29993+0.40490? 

-0.55384+0.46277? 

0.58139+0.37949? 



Teukolsky, s = —2, n = 1 



TO = -2 



m = — 1 



TO = 



m = 1 



TO = 2 



0.00 -0.081135+0.067726i -0.081135+0.067726? 
0.20 -0.042052+0.057934? -0.061317+0.063041i 
0.40 -0.023703+0.046853? -0.053312+0.056703? 
0.50 -0.018645+0.042289? -0.052550+0.053249? 
0.60 -0.015312+0.038467? -0.053615+0.049396? 
0.80 -0.012449+0.033074? -0.061191+0.038443? 
0.90 -0.012978+0.031609i -0.068071+0.029129't 
0.98 -0.015950+0.031612? -0.077430+0.015783? 



-0.081135+0.067726i 
-0.086857+0.065808i 
-0.10514+0.05697? 
-0.11959+0.04661? 
-0.13714+0.02823? 
-0.16317-0.05596i 
-0.13659-0.13081't 
-0.05976-0.18580? 



-0.081135+0.067726i 
-0.12023+0.06413i 
-0.18937+0.02718? 
-0.23424-0.02592i 
-0.27023-0.12799? 
-0.00912-0.52799? 
0.60114-0.32896't 
0.14397+0.42835? 



■0.081135+0.067726? 
-0.16306+0.05488i 
-0.30908-0.07767? 
-0.36754-0.28206? 
-0.26075-0.66420? 
2.1447-0.0490? 
-2.2114+3.6692? 
3.7992+4.8108? 
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TABLE VI: Amplitudes Aout of the first two overtones with s = —2 and I = 2, for several values of the rotation parameter and 
all values of m. 

Sasaki-Nakamura, s = —2, n = 



-1 







0.00 -4.7608+2.3771i 

0.20 -6.2957+3.76642 

0.40 -7.8741+5.8832i 

0.50 -8.5905+7.3807i 

0.60 -9.1362+9.34694 

0.80 -8.311+16.0351 

0.90 -3.399+22.656i 

0.98 23.657+28.208i 



-4.7608+2.3771't -4.7608+2.3771't -4.7608+2.3771't -4.7608+2.3771't 

-5.4365+3.0202't -4.6676+2. 3909't -3. 9841+1. 8690i -3.3819+1.4446't 

-5.9349+3.8929i -4.3699+2.4414i -3.1373+1.4490i -2.2068+0.8273i 

-6.0631+4.4740? -4.1247+2.4901i -2.6907+1.2641i -1.7053+0.6084i 

-6.0428+5.21331 -3.7910+2.5653? -2. 2271+1. 0932i -1.2657+0.4363i 

-4.8599+7.6211i -2. 6288+2. 8690i -1.2248+0.7843i -0.58015+0.18762j 

-2.0780+9.8297't -1.3270+3.1477't -0.63687+0.62642j -0.32078+0.07721i 

9.6770+9.8318? 2.1270+2.4156? 0.03281+0.34446? -0.09492-0.12174? 



Teukolsky, s 







m = -2 



m — —1 



m = 



m = 2 



0.00 
0.20 
0.40 
0.50 
0.60 
0.80 
0.90 
0.98 



0.39793 
0.38459-1 
0.35498- 

0.32915- 
0.28965- 
0.10472- 
-0.18616 
-i.2217 



0.36149? 
0.43342? 
0.51546? 

0.56460? 
0.62309i 
0.79207't 
-0.91652/" 
-0.70 Kit 



0.39793- 
0.38883- 
0.36539- 

0.34451- 
0.31268- 
0.16769- 
-0.05125 
-0.77025 



0.36149? 
0.39777? 
0.44015? 

0.46621? 
0.49797? 
0.59291? 
-0.659407 
-().i:!i9i/: 



0.39793- 
0.39300- 
0.37598- 

0.36044- 
0.33704- 
0.23613- 
0.09557- 
-().:!i9."):)- 



0.36149? 
0.36207? 
0.36462? 
0.36748? 
0.37241? 
0.39539? 
0.415187 
-().27()()7/ 



0.39793- 
0.39696- 
0.38528- 
0.37380- 
0.35645- 
0.28646- 
0.20434- 
0.02188- 



0.36149? 
0.32633? 
0.28888? 

0.26837? 
0.24626? 
0.19715i 
0.172027 
0.13092/ 



0.39793- 
0.40050- 
0.39152- 

0.38051- 
0.36232- 
0.28191- 
0.19003H 

0.01917n 



0.36149? 
0.29056? 
0.21308? 

0.16934? 
0.12071? 
0.00589? 
-0.05215? 
-().05:!37't 



Sasaki-Nakamura, s = —2, n = 1 



-1 







m ■■ 



m ■■ 



0.00 2.2945 -7.1873? 2.2945-7.1873? 2.2945-7.1873? 2.2945-7.1873? 2.2945-7.1873? 

0.20 2.5658 -9.0390? 2.3970-8.0165? 2.2247-7.0665? 2.0424-6.1893? 1.8441-5.3859? 

0.40 2.543 -10.849i 2.2599-8.6166? 1.9888-6.6818? 1.6731-5.0445? 1.2672-3.7148? 

0.50 2.344 -11.695i 2. 0427-8. 7930i 1.7817-6.3679i 1.4444-4.4156't 0.9446-2.9579? 

0.60 1.933 -12.458i 1.6641-8.8477? 1.4852-5.9476? 1. 1812-3. 7433't 0.6120-2.2745? 

0.80 -0.205 -13.360? -3.1529-8.2731? 0.4135-4.5896? 0.5162-2.2182? 0.0351-1.2045? 

0.90 -3.006 -12.697? -1.9831-7.0189? -0.6196-3.3073i 0.0900-1.2838? -0.04965-0.79665? 

0.98 -8.7106 -7.3992? -5.0880-2.5636? -1.8011-0.6338? -0.19532-0.23025? 0.07333-0.27876? 



0.00 0.82322+0.75284? 
0.20 0.84911+0.70754? 

0.40 0. 85814+0. 65292i 
0.50 0.85843+0.61654i 
0.60 0.85685+0.56852i 
0.80 0.84873+0.39392? 
0.90 0.82922+0.192491 
0.98 0.65127-0.27305? 



Teukolsky, s = — 2, n = 1 



0.82322+0.75284? 
0.77073+0.79838? 
0.67346+0.84688? 
0.59731+0.86954? 
0.49254+0.88500? 
0.14457+0.81785? 
-0.10621+0.60268? 
-0.17382+0.14948? 



0.82322+0.75284? 
0.82947+0.72453? 
0.81337+0.67702'i 
0.79766+0.64004i 
0.77713+0.58786i 
0.72054+0.39027? 
0.66918+0.17074? 
0.44118-0.27437? 



0.82322+0.75284? 
0.80991+0.74533? 
0.76818+0.71734i 
0.73512+0.68991? 
0.69265+0.64679i 
0.57393+0.46234? 
0.48666+0.24733i 
0.27060-0.14913? 



0.82322+0.75284? 
0.79038+0.76996? 
0.72213+0.77414? 
0.66983+0.76696? 
0.60121+0.74743? 
0.39631+0.61795? 
0.25229+0.42940? 
0.10617+0.07862? 



TABLE VIL Values of the QNM frequencies and of the angular momentum parameter v in the Schwarzschild limit. 
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.s = 0, / = 2 


V 

,s = -L / = 1 


.s = -2, / = 2 




1 


1.64072+0.24770? 
-1.91021-0.44483i 


1.85364+0.20371? 
-1. 17974-0. 29072i 


1.72837+0.23396? 
-2.00310-0.37593i 


n 


s = 0,1 = 2 


s = -\,l = \ 


s = -2, Z = 2 



1 


0.48364-0.09676? 
0.46385-0.29560? 


0.24826-0.09249? 
0.21452-0.29367? 


0.37367-0.08896? 
0.34671-0.27392? 



